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Abstract 



The tau lepton lifetime has been measured with the e + e — > t + t events col- 
lected by the DELPHI detector at LEP in the years 1991-1995. Three dif- 
ferent methods have been exploited, using both one-prong and three-prong r 
decay channels. Two measurements have been made using events in which 
both taus decay to a single charged particle. Combining these measure- 
ments gave r T (l prong) = 291.8 ± 2.3 s t a t ± l-5 sys fs. A third measurement 
using taus which decayed to three charged particles yielded r T (3 prong) = 
288.6 ± 2A stat ± 1.3 sys fs. These were combined with previous DELPHI re- 
sults to measure the tau lifetime, using the full LEP1 data sample, to be 
t t = 290.9 ± lA stat ± 1.0 sys fs. 
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1 Introduction 

The tau lepton is a fundamental constituent of the Standard Model and its lifetime 
can be used to test the model's predictions. In particular, lepton universality can be 
probed using the relationships 

T * = UJ BR (r —> e »*) • (D 

\g T J \m T J » f(ml/m 2 T )r T RC 

where r M)7 - and m^ jT are the lifetimes and masses of the muon and tau lepton, g e ,n, T are the 
coupling constants to the W ± for the electron, muon and tau respectively, / are phase 
space factors and r^ c are radiative corrections to the decay widths [1]. To the precision 
with which the tau lifetime and branching ratios can be measured, f(m 2 /m 2 ) = 0.9726 
while /(ml/m 2 ^) and f(m 2 /m 2 ) are 1.000; the electroweak radiative corrections r RC and 
r T RC amount to 0.9956 and 0.9960 respectively. 

The lifetime measurements presented here were performed with the data taken by 
the DELPHI experiment at the LEP electron-positron collider at centre-of-mass energies 
of the e + e~ system around 91 GeV, where tau leptons were pair-produced through the 
decay of the Z boson. As in previous measurements [2], the three-layer silicon microvertex 
detector [3] and its excellent spatial resolution were the key to achieving the precision on 
track measurements necessary to determine the short tau decay distance. 

Three techniques were used to measure the lifetime depending on the final-state topol- 
ogy of the event. In the channel in which a tau decayed into a final state containing three 
charged particles (3- prong decays), it was possible to reconstruct the decay vertex and 
measure the decay distance from the centre of the interaction region of the LEP beams. 
An analysis of the complete sample of such decays collected by DELPHI from 1991 to 
1995 was performed. When a tau decayed into final states with only one charged particle 
(one-prong decays), the lifetime information was contained in the impact parameter of 
that particle with respect to the centre of the interaction region. Two complementary 
methods, similar to the ones used on 1991-1993 data [2,4] were applied to the data col- 
lected during the 1994 and 1995 LEP running. These methods exploited the correlation 
between the impact parameters of the two charged particles in two one-prong tau decays. 

All these methods measured the tau decay length. The conversion to a lifetime 
used the Lorentz boost parameter 7/5 which was estimated from the tau mass m T = 
1776.99±g;|§ MeV/c 2 [5] and the energy of the LEP beams. 

The Monte Carlo program KORALZ 4.0 [6], together with the TAUOLA 2.5 [7] li- 
brary were used to model tau-pair production and decay. Backgrounds were studied 
using several generators: DYMU3 [8] for e + e~ — * events; BABAMC [9] and BH- 

WIDE [10] for e + e _ — > e + e~ events; JETSET 7.3 with specially tuned fragmentation 
parameters [11] for e + e~ — > qq events; BDK [12] for reactions with four leptons in the 
final state, including two-photon events where one or two e + or e~ were not observed in 
the detector; and TWOGAM [13] for e + e~ — > (e + e~)qq events. The generated events 
were interfaced to a detailed model of the detector response [14] and reconstructed with 
the same program as the real data. Separate samples were produced corresponding to 
the detector configurations in different years. 

The DELPHI detector is described in [15]. This analysis used the charged particle 
tracking system covering the polar angle range |cos#| < 0.73. This consisted of four 
detectors in a 1.2 T solenoidal magnetic field. 
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• The Microvertex Detector (VD) was a three-layer silicon vertex detector, which 
provided an Rcf) 1 precision of 7.6 /im and a two-track separation of 100 /im. In 
1994 an upgraded version with two of the three layers equipped with double-sided 
detectors was installed, providing a z precision of 9 for tracks perpendicular to 
the beam direction. 

• The Inner Detector (ID) was a gas detector with a jet-chamber geometry. It mea- 
sured up to 24 Bxp coordinates per track, yielding a track element with an B(f) 
precision of 50 /im. 

• The Time Projection Chamber (TPC) was the main tracking detector of DELPHI, 
situated between radii of 30 cm and 120 cm. Up to 16 points per track produced a 
track element with an Bxp precision of 250 /im. 

• The Outer Detector (OD) consisted of 24 modules containing 5 layers of drift tubes 
operating in limited streamer mode and situated at a radius of 2 m. Charged particles 
produced track elements with 300 /im precision in Bxp. 

The most important figure of merit in the tau lifetime measurement is the precision on 
the impact parameter defined as the distance of closest approach of a track extrapolated 
to the centre of the interaction region. The impact parameter was given the same sign as 

— * — * 

the z component of d x pr where d is the projection on the B(f) plane of the vector from 
the centre of the interaction region to the point of closest approach and pr the projection 
on the same plane of the particle momentum. 

The precision on the impact parameter was extensively analyzed in the detailed study 
of tracking uncertainties for the DELPHI measurement of the Z decay width into pairs 
of bb quarks [16]. Typically these errors correspond to an uncertainty on the impact 
parameter of 

65 /xm 

^track = 20 /im © . (3) 

p[GeV/c] sm J/ ^ 9 

where © indicates the sum in quadrature. 

Nevertheless this estimation was not fully satisfactory to describe the impact parame- 
ter resolution in tau decays for two reasons. Firstly, the topology of b events differed from 
that of tau decays. Ambiguities in the reconstruction of overlapping tracks degraded the 
resolution in b events. This did not affect one-prong tau decays but was an even more 
severe problem in three-prong decays because of the small opening angle of the tau decay 
products. Secondly, particles from tau decays contain a considerably higher fraction of 
muons and electrons than particles from b-decays. This affected the tracking precision in 
different ways since muons had no hadronic scattering in the detector material leading to 
a better tracking precision, while electrons were strongly affected by the bremsstrahlung 
process leading to a worse precision. Therefore the analyses needed to perform checks of 
the precision in the specific topology under study. 

Since it was not possible to determine the production point of the tau pair, the centre 
of the interaction region was used. An error was induced by this approximation due to 
both the size of the beams and the accuracy of the estimated position of the centre of the 
interaction region, which was reconstructed from the distribution of the primary vertices 
of hadronic Z decays. In the vertical (y) direction the beams were very narrow and the 
main contribution came from the uncertainty on the vertex reconstruction, which resulted 
in a typical precision of 10 /im in the production point position. In the horizontal (x) 
direction, the uncertainty was dominated by the beam size and depended on the LEP 

1 R, <j> and z define a cylindrical coordinate system, +z being coincident with the electron beam and _R and (p in the 
plane transverse to the beam. The angle 8 is the polar angle defined with respect to the z axis. 
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operating conditions. The resulting uncertainty on the production point position ranged 
from about 90 /im in 1992 to 160 /im in 1995. 

In addition to the above mentioned tracking detectors, the identification of the r decay 
products relied on electromagnetic calorimetry for electron identification, and on hadron 
calorimetry and muon chambers for muon identification. 

• The barrel electromagnetic calorimeter was a High-density Projection Chamber 
(HPC), covering the polar angle region from 43° to 137°. It had a high granularity 
and provided nine layers of sampling of shower energies. 

• The Hadron Calorimeter (HCAL) was situated outside the magnet solenoid and had 
a depth of 110 cm of iron. It was sensitive to hadronic showers and minimum ionizing 
particles and consisted of four layers with a granularity of 3.75° in polar angle and 
2.96° in azimuthal angle. 

• The barrel Muon Chambers (MUB) consisted of two layers of drift chambers, the first 
one situated outside 90 cm of iron and the second outside the hadron calorimeter. 
The acceptance in polar angle of the outer layer was slightly smaller than the other 
barrel detectors and covered the range | cos#| < 0.602. The polar angle range 0.743 < 
| cos#| < 0.940 was covered by the forward Muon Chambers (MUF). In 1994 a layer 
of Surround Muon Chambers (SMC), based on limited streamer tubes, was installed 
to fill the gap between the barrel and forward regions. 

• The TPC also provided up to 192 ionisation measurements per charged particle 
track, which were useful for electron/hadron separation. 

Section 2 describes the decay length analysis applied to three-prong tau decays. The 
two one-prong measurements, together with the correlations induced from using overlap- 
ping data samples, are described in Section 3. Section 4 presents the combined result for 
the full 1991-1995 data sample. Conclusions are reported in Section 5. 

2 The Decay Vertex Method 

This method is an improved version of the one described in [2]. Decays of the Z into 
two taus were selected, where both taus decayed to three charged particles plus neutral 
particles (the 3v3 topology), or where one tau decay contained three charged particles 
and the other decay one charged particle (the 3vl topology). In the 3vl topology only 
the three-prong decays were used. 

The event was divided into two hemispheres defined by the plane perpendicular to 
the thrust axis and passing through the centre of the interaction region. The highest 
momentum charged particles in the two hemispheres were required to be separated by an 
angle of at least 2.9 and 2.975 radians in the 3vl and 3v3 topologies respectively. The 
tighter cut for the 3v3 topology reflects the increased background from hadronic decays 
of the Z. For each tau decay, the three charged particle tracks had to be accurately 
reconstructed: each track had to be associated to at least two microvertex detector hits 
and all three tracks had to be consistent with coming from a common vertex with a \ 2 
probability greater than 0.5%. If the momentum of any charged particle in an event was 
greater than 35 GeV/c, there had to be no indication that this particle was either an 
electron or a muon. This was necessary in order to eliminate four-fermion events with 
e + e~e + e~ and e + e~/z + /i~ final states. After this selection the only appreciable background 
consisted of e + e~T + T~ events and hadronic Z decays. 

From simulation the background fractions, 6j vl and 6j v3 , of four-fermion e + e _ r + T~ 
events remaining in the 3vl and 3v3 topologies were estimated to be (0.15 ± 0.03)% 
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and (0.5 ± 0.2)% respectively. These events had a decay length which was on average 
25% smaller than those produced directly, due to the lower effective r + r~ centre-of-mass 
energy. This corresponded to a bias of —0.1 fs in the determination of the lifetime. 

The background from hadronic Z decays was reduced by requiring at most six elec- 
tromagnetic neutral deposits in the event. By relaxing the cut on the opening angle and 
comparing the amount of additional events in data and simulation, the hadronic contam- 
inations in the 3vl and 3v3 topologies, b 3 ? 1 and b 3 ^ 3 were estimated to be (0.38 ± 0.06)% 
and (0.8 ± 0.2)% respectively. These contributions were included in the lifetime fit with 
their errors taken as systematic uncertainties. The hadronic background had no life- 
time content and corresponded to a bias of —1.2 ± 0.2 fs (—2.2 ± 0.5 fs) in the lifetime 
determination in the 3vl (3v3) topology. 

Within the 1991-1995 data sample, 15427 and 2101 three-prong tau decays were se- 
lected in the 3vl and 3v3 topologies respectively. The distance in the transverse plane 
between the reconstructed decay point of the tau and its production point, taken as the 
centre of the beam crossing region, was calculated taking into account the reconstruction 
uncertainty of the former and the size of the latter. This was converted to a decay time 
t and associated uncertainty a by dividing by ^f3(^fs)c sin 9 where 9 is the polar angle 
of the thrust axis of the event, c is the speed of light and 7/3(-\/s) — \/s/Arn* — 1, with 
\fs the centre-of-mass energy of the collision, taken as twice the beam energy with a 
correction for initial and final state radiation that was evaluated from the simulation to 
be between 0.8% and 1.4% depending on the beam energy. The accuracy in the energy 
correction resulted in a systematic uncertainty of 0.1 fs on the lifetime determination. 

Simulation studies showed that the central value of the pull distribution (defined as 
the reconstructed decay time less the true decay time, divided by the uncertainty) had 
a slight positive shift of 0.02 ± 0.01, leading to a lifetime bias of 1.6 ± 0.8 fs. This bias 
was subtracted from each event with the uncertainty taken as a systematic uncertainty 
on the final lifetime determination. 

The tau lifetime, r r , was extracted from a maximum log-likelihood fit to the data 
using the function 

L(t t , A£ r , Af, A yr ) = log P {U\a h r T , A£ r , Af , \%) , yr=1991-1995, (4) 

where ti is the measured proper time for event i, after applying the correction for the 
biases resulting from reconstruction and alignment of the Microvertex Detector, Oi its 
uncertainty as computed by the reconstruction program and A yr are scaling factors as 
explained below. 

The probability density function P is given by the convolution of the probability 
density function of the positions of decay vertices and the resolution function of the 
detector. The sum runs over the full sample of events, spanning the five data-taking 
periods used in the analysis. 

The probability density function of the decay vertices contains the dependence upon 
the lifetime: 

/vtx(*|r r ) = (1 - b)° p ° - C° WK) + &y po £(t|0.75r r ) + C°<H0, topo=3vl, 3v3, (5) 

where E(x\r]) denotes a normalized exponential distribution with a decay constant rj and 
5(t) is a Dirac delta function. The factor 0.75 in the four-fermion background term takes 
into account the reduced centre-of-mass energy for these events. 

The lifetime component was convoluted with a resolution function / ros which, accord- 
ing to the simulation, could be adequately parameterized as the sum of three Gaussian 
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Year 


Scale factors 




Ao 


Ai 


A 2 


1991 


1.13 ±0.29 


1.30 ±0.25 


1.66 ±0.62 


1992 


0.79 ±0.51 


0.87 ±0.21 


0.86 ±0.17 


1993 


1.03 ±0.18 


0.86 ±0.10 


1.08 ±0.17 


1994 


1.15 ±0.07 


0.99 ±0.05 


1.43 ±0.18 


1995 


1.06 ±0.17 


0.93 ±0.10 


1.21 ±0.22 



Table 1: Fitted values for the resolution function parameters in (7). 



distributions: 

fr es (t\a, / 2 , fs, h, k 2 , k 3 ) = (1 - h ~ fs)G(t\ha) + f 2 G(t\k 2 a) + fcG(t\k 3 a), (6) 

where G{x\if) is a normalized Gaussian function centred at zero and with width rj. f 2 
and f 3 are the fractions of the second and third Gaussian functions, and k\a, k 2 a, k 3 a are 
the widths. The shape of the resolution function was taken from the simulation, where 
the three Gaussian functions had widths of k\ = 0.97, k 2 = 1.6, k 3 = 5.1 times a and the 
fractional contributions of the second and third Gaussians were f 2 = 0.25 and f 3 = 0.007 
respectively. However, there was some indication that the widths and proportions of 
the Gaussian functions were slightly different in the real data and, moreover, that they 
varied from year to year. Three scale factors for each year of data-taking, Ag r , Af , \ y 2 r 
were introduced in the likelihood function and these affect the fraction, f 2 , of the second 
Gaussian and the widths of the first and second Gaussian functions, ki,k 2 . 
The complete expression for P is given by 

P (t\a u r T , Xl\ Af , A^ r ) = Mt\r T ) ® / res (*K (A^), fs, (X?h), (A^ 2 ), fc 3 ). (7) 

The scale factors were fitted together with the lifetime and their values are reported in 
Table 1. The tau lifetime was estimated to be 288.6 ±2.4 fs, and most of the scale factors 
were seen to be consistent with unity, confirming that the shape and parameterization 
of the resolution function in the simulation were reasonable. The observed decay length 
distribution compared to the result of the maximum likelihood fit is shown in Fig. 1. 

A number of consistency checks were performed. Firstly, the maximum likelihood fit 
was repeated holding all the scale factors fixed to unity (i.e. taking the resolution function 
from the simulation) and this gave a value of r T = 289.1 fs. Secondly, the resolution 
function was taken from the simulation with a single scale factor allowed to multiply 
ki,k 2 and k 3 . This gave a value of r r = 289.0 fs. Thirdly, starting with the simulation 
resolution function, five scale factors (one for each year) were allowed to multiply f 3 , in 
order to gauge the effect of the smallest but broadest Gaussian. This gave a value of 
r T = 288.0 fs. Fourthly, a weighted mean 2 of the data, corrected for background biases 
was computed. This test is insensitive to the exact shape of the distribution and it gave 
a value of r T = 288.0 fs. 

The largest systematic uncertainty came from the accuracy of the alignment of the 
Microvertex Detector. This was calculated year-by-year by taking a sample of hadronic 
decays of the Z with three tracks in one hemisphere and more than three tracks in the 
other. In such topologies, the momentum and invariant mass distributions of the three 

2 The weighted mean is defined as Ei(W^<)/ E { l / w i) where Wi = (<rf MS ) 2 + r 2 and af MS is the RMS of the 
resolution function for a given Uj. 
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Year 


Alignment correction 


Contr. to syst. uncertainty 


1991 


83 ± 45 fan 


±0.23 fs 


1992 


-2 ± 19 jum 


±0.44 fs 


1993 


0± 17/im 


±0.42 fs 


1994 


21 ± 13 /im 


±0.57 fs 


1995 


32 ± 22 /im 


±0.47 fs 


Total 




±0.98 fs 



Table 2: Alignment corrections from hadronic data and their contributions to the total 
systematic uncertainty. 

tracks were similar to those in tau decays. A vertex was formed from the three tau- 
like tracks and the distribution for the decay distance was compared to a simulation of 
hadronic decays of the Z. Possible indications for shifts in the reconstructed position 
were observed and are summarized in Table 2. These corrections were applied to the 
reconstructed decay distances before the maximum likelihood fit was made. The related 
systematic uncertainty was calculated for each year by changing the alignment correction 
for that year by its uncertainty and repeating the complete fit using all years. The align- 
ment corrections were uncorrelated between years and so their systematic uncertainties 
were added in quadrature to give ±0.98 fs. 

A summary of the systematic uncertainties due to each source is given in Table 3. The 
tau lifetime from three-prong decays of the tau lepton was thus measured to be 

r r (3 prong) = 288.6 ± 2A stat ± l.3 sys fs. 

This result supersedes previously published three-prong DELPHI results. 

3 One-Prong Lifetime Measurements 

The lifetime information from one-prong tau decays was obtained by measuring the 
impact parameters of the charged decay products. In the case of perfect knowledge of 
the track parameters and of the production point, the impact parameter in the R<p plane 
is given by: 

d = Lsin6> T sin(0 — <p T ), (8) 

where L is the decay length, <p T the azimuthal direction of the decaying object, the 
track's azimuth and 6 T the polar angle of the decaying object. Impact parameters, signed 
according to the geometrical definition given in Section 1, were used in the calculation of 
the resolution functions, as well as in the extraction of the tau lifetime. 

The tau lifetime was extracted from events in which both taus decayed to a sin- 
gle charged particle using two methods. The first method used the impact parameter 
difference (IPD) which represents an improvement over the single hemisphere impact 
parameter lifetime determination (see for example [4]) by reducing the dependence of 
the measurement on the unknown tau-decay angle. In the IPD method, knowledge of 
the tau-pair production point was limited by the size of the interaction region, whose 
dimensions are larger than the track extrapolation precision of the detector. To avoid 
this limitation, the second method used the track pair miss distance (MD). In the MD 
method the two impact parameters in a t + t~ event were summed so that the dependence 
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Fitted lifetime 


288.6 ± 2.4 fs 


Background 


±0.2 fs 


Radiative Energy Loss 


±0.1 fs 


Reconstruction Bias 


±0.8 fs 


Alignment 


±1.0 fs 


Total 


288.6 ±2.4 ±1.3 fs 



Table 3: Summary of fit result and systematic uncertainties for the 3-prong topology. 

on the production point inside the interaction region cancelled to first order. This second 
method was sensitive to the knowledge of the resolution function, which is described in 
detail in Section 3.3. 

The t + t~ data were selected in the same way as those used for the tau polarisation 
measurement [17]. The event was divided into two hemispheres by a plane perpendicular 
to the thrust axis and events were required to have: 

• the highest momentum charged particle in at least one of the two hemispheres with 
|cos#| < 0.732; 

• total charged particle multiplicity less than 6; 

• isolation angle between the highest momentum particles in the two hemispheres 
greater than 160°; 

• the highest momentum particles in the two hemispheres passing at less than 4.5 cm 
in z and less than 1.5 cm in the R<p plane from the center of the interaction region 
and the difference in z of the points of closest approach should be less than 3 cm; 

• total energy of the event greater than 8 GeV and total transverse momentum greater 
than 0.4 GeV/c; 

• acollinearity of the two highest momentum tracks in each hemisphere greater than 
0.5°; 

• P rad less than the beam momentum and P rad less than the beam energy. 

In the above P rad = (\pi\ 2 + Ij^I 2 ) 1 ^ 2 and p\ and P2 are the charged-particle momenta 
and P rad = (Ef + El) 1 / 2 . The variables E\ and E 2 are the total electromagnetic energies 
deposited in cones of half-angle 30° about the charged-particle momentum vectors p\ and 
P2 respectively. 

In addition there must be only one charged particle track per hemisphere with hits 
in the Microvertex Detector. Both tracks had to satisfy the following requirements: a 
momentum transverse to the beam axis greater than 1 GeV/c, associated R<p hits in at 
least two VD layers, at least 11 associated points in the TPC and a x 2 probability about 
0.1% for the additional x 2 when all of the VD hits are added to the track fitted to the 
TPC segments. 

After the application of the above criteria, the remaining background came essentially 
from e + e~, n + fi~ and 77 — > il events. Further selections were applied depending on 
the result of the lepton identification described in Section 3.1. To suppress fi + fi~ events, 
if both charged particles were identified as muons, P rad < 35 GeV/c was required. To 
suppress e + e~ events, if one of the particles was identified as an electron and the other 
was not a muon, the requirements P rad < 35 GeV/c and P rad < 30 GeV were imposed. 
Finally two-photon events were suppressed requiring P rad > 11 GeV/c or E rad > 8 GeV, 
if both charged particles were identified as muons or electrons. 
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These criteria selected 17366 and 8670 events in the 1994 and 1995 data samples respec- 
tively. The residual background in the sample was obtained from simulated background 
events and amounted to (0.58 ± 0.05)% from dilepton events and to (0.31 ± 0.03)% and 
(0.46 ± 0.04)% in 1994 and 1995 respectively from 77 events. The quoted uncertainties 
were derived from the available statistics of simulated events. 

3.1 Lepton Identification 

Particle identification had an influence on different aspects of this analysis, namely 
background suppression and, for the impact parameter sum method, also in the physics 
and resolution function parameterization (see Section 3.3). A neural network was imple- 
mented to improve the particle identification, using the same variables as for the particle 
identification in r decay applied in [18]: 

• the neutral electromagnetic energy measured by the HPC in a cone of half-angle 19° 
about the track; 

• the number of hits associated in the first HCAL layer; 

• the number of hits associated in the last HCAL layer; 

• the average energy deposited in a HCAL layer; 

• the number of hits associated in the muon chambers; 

• the pull functions of the measured energy loss dE/dx compared with the one ex- 
pected from an electron or from a pion 3 ; 

• the E/p of the particle, where p is the momentum measured from the curvature of 
the track and E is its associated electromagnetic energy; 

• the x 2 probability for dz = z H pc — z extI , where Zhpc is the z of the shower associated 
to the track, measured by the HPC and z cxtv the one extrapolated from the TPC to 
the HPC. 

A different neural network was implemented for the small fraction (2.5%) of events with no 
reliable dE/dx measurements. Both neural networks were three-layer feed-forward neural 
networks, with only two output nodes, one for electrons and one for muons. Hadrons were 
classified as all charged particles not identified as either electron or muon. 

The neural network performance was tested with simulated samples different from 
those used for the training. For a simulated sample corresponding to 1995 data, an 
electron identification efficiency and purity of (95.6±0.1)% and (91.8±0.1)% respectively 
were obtained, while for muons the corresponding numbers were (96.7±0.1)% and (95. 0± 
0.1)%. 

3.2 The Impact Parameter Difference Method 

At LEP, the taus were produced with high boost and in back-to-back pairs which 
allowed some simplifications [19] in evaluating equation (8). The boost allows the sine of 
the decay angle to be approximated by sin(</> — <p T ) m <p — T , while the collinearity allows 
the substitutions sin# T - = sin# T + and <p T - = <p T + ±tt. Therefore, considering both tracks 
in the event, 

d + — gL = sin# T + (L + (j) + — L-4>- — (L + — LJ)4> T + ± L_7r) (9) 

where the variables are the same as in equation (8), with the additional subscript to 
indicate the charge of the tau in the pair. The sign of n is chosen to normalize the 

3 The pull function is )mcas — ( ) e /,r ^ w j lere {dE / dx)meas is the energy loss measured by the TPC, (dE / 'dx) e i w 

° dE I dx ' 

is the one expected from an electron (or a pion) and a,iE /dx ls the energy loss reconstruction uncertainty. 
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acoplanarity of the two observed tracks, A0 = + — 0_ ±7r, to be in the ±7r/2 range. An 
average over the decay lengths shows that the average difference of impact parameters is 
related to the projected acoplanarity, sin# r A0, of the decay products [19]: 

(d+ - <Z_) = ~fP(Vs)cr T sin # r A0. (10) 

Equation (10) shows that the lifetime information can be extracted from the mea- 
surement of the directions of the outcoming decay products, independently of the poorly 
known production angle <f> T . The uncertainty on the impact parameter difference is, how- 
ever, affected strongly by the uncertainty on the production point. Moreover the small 
angle and tau-pair collinearity are just approximations that cease to be valid for sub- 
stantial decay angles and for hard initial-state radiation. Therefore deviations from the 
linear behaviour of equation (10) are expected at large values of acoplanarity. 

As discussed in a previous publication [2], the average decay length (L) = r yj3(Mz)cT T 
was determined from an event weighted x 2 fit °f a straight line to the impact parameter 
difference as a function of the projected acoplanarity: 

= E w E WX V ~Y. WX Y. W V / n x 

E w E wx<1 ~~ (E wx) 2 ' 

where the sum extends over the selected events, 

V = d+-d-, (12) 

X = {^) Sm9 ^ + '^- (13) 
The weight w given to each single measurement depends upon several variables: 

w = Khy S + ^x,bcam(sin0+ - sin0_) 2 + <j 2 beam (cos 0+ - cos0_) 2 + a 2 d+ + a 2 d _y l , (14) 
where 

• Ophys is the r.m.s. of the distribution of the impact parameter difference due to the 
lifetime spread; 

• °x.bcam and o"y ib eam axe the uncertainties in the production point due to the beam 
size and the knowledge of beam position in the x and y DELPHI coordinates (see 
Section 1); 

• (Td± is the impact parameter uncertainty from the DELPHI track fit. 

An outlier rejection was performed discarding the events with the highest significance 
of the residual \y/w(y — (L)x)\. This procedure, designed to reduce the statistical fluc- 
tuations in the final result due to badly reconstructed events, was used to reject at most 
1% of the sample. 

The impact parameter difference versus the projected acoplanarity is shown in Fig. 2, 
together with the result of the fit. The results of the fits to 1994 and 1995 data were: 

(L) 94 = 2161 ±33 jim, 
(L) 95 = 2150 ±51 Aim, 

where only the statistical uncertainties are shown. The x 2 /DOF Q f the fits were respec- 
tively 16279/16288 and 7544/8098 for 1994 and 1995 data, the ratio being lower than 
unity as expected due to the outlier rejection. 

The IPD method suffered from several biases which had their origin in the assumptions 
and approximations of the method and in the need to trim the tails of the significance 
distribution. Another source of bias was the background contamination. 
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The assumption that collinear tau pairs had been produced with the full centre-of- 
mass energy led to a bias of —26.9 ± 1.3 /im which was evaluated from the simulation of 
initial-state radiation. 

The above mentioned outlier rejection, performed by trimming of the residuals, in- 
troduced an additional bias since the asymmetric exponential tail due to the lifetime 
distribution was preferentially cut. The effect was significant and the induced shift was 
derived from simulated event studies. A check was made comparing the behaviour of 
the decay length fitted on the data with the expectation from the simulation. The fitted 
decay lengths as a function of the trim fraction are displayed in Fig. 3. There was good 
agreement between data and simulation, providing a lifetime determination that is stable 
with respect to changes in the trim fraction. A trimming point of 0.4% was chosen as 
in previous publications, corresponding to a bias of —43.4 ± 7.9 /jm and —40.3 ± 6.8 
/im for 1994 and 1995 respectively. The uncertainty quoted is the one due to simulation 
statistics. In order to take into account the uncertainty due to the modeling of the trim 
dependence in the simulation, a further contribution to the systematic uncertainty was 
added. It was estimated as the maximum difference between the lifetime evaluated at the 
chosen trim point and all the other evaluations in the interval [0.1;1.0]%. This amounted 
to 9 /im and 10 /im for 1994 and 1995 data respectively. 

As an additional systematic check, the point at which the line intercepted the y-axis 
was determined. This was expected to be different from zero, due to the correlation 
between the energy loss (particularly relevant for electrons) and the shift in the recon- 
structed impact parameter. This effect was checked using the lepton identification algo- 
rithm described in Section 3.1. Combining the 1994 and 1995 samples, the offsets were 
23.0 ±2.7 /im for events containing at least one identified electron and 3.9 ±1.4 /im for 
events with no identified electrons. These results agreed well with the expectations from 
the simulation of 20.6 ± 0.7 /im and 5.4 ± 0.6 /im respectively. 

The bias induced by the background was estimated by adding to a simulated sample 
50 samples of simulated background events which had passed the selection criteria. The 
average bias was —11.4 ± 4.9 /im for 1994 and —17.5 ± 7.4 /im for 1995, where the 
systematic uncertainty was the r.m.s. of the biases calculated in each of the 50 samples. 
Additional systematic uncertainties were due to the uncertainty on the resolution function 
(3.8 /im) and to the vertex detector alignment. The latter was checked by computing 
the lifetime using a vertex detector geometrical description simulating the alignment 
uncertainties. Most parameters are well constrained by the alignement procedure and 
provide negligible variation in the lifetime determination. Only the less constrained 
deformation, a coherent radial variation of the silicon ladders, provided a visible shift 
in the reconstructed values, which amounted to 3.1 /im, for a radial movement of 20 
/im. These results are in qualitative agreement with what obtained from a simplified 
simulation in section 5.3 of [20]. 

The uncertainty due to the resolution function was considered to be correlated between 
the two years and a correlation of 0.45 was calculated for the uncertainty due to the 
background. All the other systematic uncertainties were uncorrected, as they came 
from calibrations which were computed with different data sets for the two samples. A 
summary of all systematic uncertainties is given in Table 4. The lifetimes were measured 
to be: 

r r (IPD, 94) = 291.7 ± 4.3 stat ± 1.8 sys fs, 
r r (IPD, 95) = 290.7 ± 6.7 stat ± 2.0 sys fs, 
r T (IPD, 94±95) = 291.4 ± 3.6 stat ± 1.5 sys fs. 
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1994 


1995 


Correlation 




[fs] 


r r 1 

[fs] 


coefficient 


Fit result 


281.1 ± 4.3 


279.7 ± 6.7 





byst. sources: 








TV yr l1 l l * 

Method bias 


i o r i no 

+3.5 ± 0.2 


+3.5 ± 0.2 


1 


trim 


+5.6 ± 1.0 


+5.2 + 0.9 





Trim data/MC agreement 


±1.2 


±1.3 





Background 


+1.5 + 0.6 


+2.3 + 1.0 


0.45 


Alignment 


±0.4 


±0.4 





Resolution 


±0.5 


±0.5 


1 


Result 


291.7 ± 4.3 ± 1.8 


290.7 ±6.7 ±2.0 


0.02 


Average 94+95 


291.4 ± 3.6 ± 1.5 



Table 4: Summary of results, correction for biases, systematic uncertainties and combi- 
nation of the Impact Parameter Difference measurements. 



3.3 The Miss Distance Method 

At LEP, the algebraic sum of the impact parameters, 5 = d + + cL, named the "miss 
distance", was strongly correlated to the separation of the two tracks at the production 
point. The width of the miss distance distribution for one- prong versus one-prong tau 
decays depends on the value of the lifetime. This was measured by an unbinned maximum 
likelihood fit to the observed distribution. The probability density function was given by 
the convolution of a physics function and a resolution function. 

The physics function was given by the distribution of miss distances expected from 
the decay length of r's. This was built from the convolution of the impact parameter 
distribution of tracks originating from r + decays with that of tracks originating from r _ 
decays. 

The dependence of the impact parameter distribution upon the track momentum, the 
decay kinematics and the r helicity was also considered. The differences in the distri- 
bution for leptons and hadrons depend on the different r decay kinematics in one-prong 
topologies; while leptonic tau decays have three final-state particles, hadronic decays have 
two or more depending on the number of neutral particles present in the decay. These 
differences are important for momenta less than 20 GeV/c. 

The shape of the single impact parameter distribution, as a function of decay dynamics 
and kinematics, was determined on a sample of 185842 simulated events which passed 
all the selection and quality cuts; it was parameterized as a linear combination of three 
exponentials. 

To obtain the miss distance, the convolution of the two functions describing the single 
impact parameter distribution, fjp + (d + ) for the r + and //p_(cL) for the r~ , was calcu- 
lated. Since the helicities were not known for a single event, the function was calculated 
as the sum of the physics function for positive helicity events and the one for negative 
helicity events, mixed according to the r polarisation [21]. 

To compute the lifetime, the physics function was convoluted with the experimental 
uncertainties on the reconstruction of the miss distance. As the lifetime information was 
only in the width of the miss distance distribution, and not in its average value, the 
increase in width induced by the reconstruction uncertainties had to be evaluated with 
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high precision. In particular it was essential to have a good modeling of the tails due to 
scattering of the particles through the apparatus. 

The resolution was determined for the single impact parameter and then convoluted 
to provide the miss distance. Hadrons and leptons required different resolution functions, 
as can be seen in Fig. 4, where the variance of the pull distribution of the reconstructed 
impact parameter is plotted against p noim , for the simulated tau sample. The pull is 
defined as the difference between the generated and reconstructed impact parameter 
divided by the tracking uncertainty given by the reconstruction program. The variable 
p noim is defined as 

p ( p \ 

Pnorm = p I 2 - j , (15) 

^bcam \ ^beam / 

where p is the particle momentum. This variable was chosen because its distribution is 
almost flat for tau decays. 

Fig. 4 shows that the average value of the variance of the pull is different from unity. 
This was due to the presence of tails in the impact parameter distribution. At low 
momenta all types of particles had similar precision since the dominant effect was multiple 
Coulomb scattering while at high momenta the different interactions of the particles show 
up in a difference in resolution. The momentum dependence of the variance of the impact 
parameter pull for hadrons was much less than that for muons and electrons. This was a 
consequence of the tuning procedure used to obtain the tracking uncertainties, which was 
based largely on (both simulated and real data) samples of hadrons in multi-hadronic Z 
decays. 

To have an acceptable description of the resolution function, a three-Gaussian param- 
eterization, 

3 

frMP, 0) = /*(P» 6 ) G ( d \^ *)) > ( 16 ) 

1=1 

was used. Both the relative fractions fi and widths Oi of the Gaussians were dependent 
on the momentum and polar-angle, as explained below. 

The resolution function was calibrated on the data. This required event samples with 
a topology similar to that of one-prong tau decays but with no lifetime effect. Applying 
the same track quality cuts and lepton identification as for the tau pair selection, but 
appropriately choosing the P rad and E rad cuts, it was possible to select samples with e + e~ 
and pairs, produced by e + e~ annihilation at high energy, and by 77 interactions at 

low energy. The obtained purities where 98.2% for the high energy e + e~ pairs, 99.4% for 
the high energy pairs and 94% for the low energy lepton pairs. The approach chosen 

was to use the data for the estimation of the resolution at high and low momenta and the 
simulation to interpolate in the intermediate momentum region. The parameterization 
used was 

fi = (1 - a)fmhfth + afmifth o"i = [(1 - a)k lh + aku]a , 

f 2 = (1 - a )(l - fmhjfth + "(1 - fmi)ftu <?2 = [(1 - ac)k 2h + ak 2 i]cr Q , (17) 

h = (l-a)(l-/ th ) + a(l-/«), <7 3 = [(1-P)k 3h + Pk m ]a 0: 



where 



a 



1 - 



p sin 3/2 9 


(■ 


p sin 3/2 6 


■^beam 




-^beam 


p sin 3/2 9 




p sin 3/2 9 


-^bcam 







(18) 
(19) 
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and <r is the uncertainty estimation from the b-tagging package [16]. The parameters 
fmh, fth, hh, hh, hh and f mh f u , hi, hi, hi were determined from the miss distance 
distribution at high and low momentum respectively. As an example the fitted distribu- 
tions for the 1994 data sample are shown in Figure 5. The parameters a and (3 range 
from 1 at low momentum to at high momentum. The steepness of the change was 
controlled by the parameters a and b, derived from the simulation. For hadrons an expo- 
nential contribution was added in order to take into account the effect of elastic hadronic 
interactions. 

The lifetime was determined by an unbinned maximum likelihood fit to the observed 
distribution, where the probability density function was given by the convolution of the 
physics function and the resolution function described above. The entire procedure was 
tested on simulated samples from which a bias of —0.2 ± 0.9 fs was measured. 

Fig. 6 shows the joint distribution of the miss distance for 1994 and 1995 data, with 
the best fit superimposed. The measured lifetimes, including all the corrections, are 



The sources of systematic uncertainties are listed in Table 5. The event selection cri- 
teria were varied inducing a lifetime change of 1.1 fs in 1994 and 1.0 fs in 1995. The 
influence of the physics function and resolution function was checked by varying by ±1(7 
the parameters of the functions parameterization, taking into account correlations. A 
further contribution to the resolution function, due to hadronic scattering, was evaluated 
comparing the data and simulation for hadronic events. The residual lepton misidentifi- 
cation, after applying the procedure in section 3.1 resulted in a systematic uncertainty 
of ±0.2 fs. The effect of background from e + e~, and 77 events was evaluated 

using simulated samples that passed all the selection cuts, resulting in estimated biases 
of -0.6 ± 0.4 fs in 1994 and -0.8 ± 0.4 fs in 1995. The contribution to the systematic 
uncertainty due to the alignment of the Microvertex Detector, estimated using the same 
procedure as in section 3.2, was 0.5 fs. The dependence on the tau polarization P T was 
checked by varying it in the range [—0.11; —0.17] resulting in a ±0.1 fs variation on the 
lifetime, while the effect of the transverse polarisation correlation was evaluated as ±0.4 
fs [2]. The fit was performed over the range |<5| < 1.5 mm. This range was varied between 
\5\ < 1 mm and |<5| < 2 mm to study the systematic effect. The maximum difference 
observed was taken clS du systematic uncertainty. 

The measurements for the two years were combined, accounting for the correlations 
in the systematic uncertainties shown in Table 5, to give the result 



4 Combination of Measurements 

The lifetime of the tau has been measured with three methods. The two one-prong 
measurements were performed on the same data sample and were combined taking into 
account correlated statistical and systematic uncertainties. The statistical correlation 
was obtained by subdividing the simulated data into 89 samples of 5000 events, applying 
the two analysis methods on each sample and computing the correlation coefficient as: 



r r (MD, 94) = 292.5 ± 2.8 stat ± 2.3 sys fs, 
r r (MD, 95) = 291.0 ± 4.(W ± 2.3 sys fs. 



r(MD, 94-95) = 292.0 ± 2.3 stat ± 2.1 sys fs. 



P 



(20) 
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1994 


1995 


Correlation 




[fs] 


[fsl 

L J 


coefficient 


Fit result 


291.7 ± 2.8 


290.0 ±4.0 


0.0 


Svst Sources - 








Method bias 


+0 2 ± 9 


+0 2 ± 9 


1.0 


Event Selection 


±1.1 


±1.0 


0.0 


Physics Function 


±0.8 


±0.8 


1.0 


Resolution Function 


±1.3 


±1.5 


0.9 


Particle Misidentification 


±0.2 


±0.2 


1.0 


Background 


±0.6 ±0.4 


±0.8 ±0.4 


0.7 


Alignment 


±0.5 


±0.5 


0.0 


Polarisation 


±0.4 


±0.4 


1.0 


Fit Range 


±0.7 


±0.2 


0.0 


Result 


292.5 ±2.8 ±2.3 


291.0 ±4.0 ±2.3 


0.17 


Average 94+95 


292.0 ±2.3 ±2.1 



Table 5: Summary of results, correction for biases, systematic uncertainties and combi- 
nation of the Miss Distance measurements. 



where r IPD j and r MD j are the determined lifetimes for each sample respectively with the 
impact parameter difference and the miss distance methods, and r T is the simulated 
lifetime. The resulting statistical correlation was 36%. Among the systematic uncertain- 
ties only the background estimation and alignment contribution were correlated. This 
provides a combined result of 

r T (l-prong, 94±95) = 291.8 ± 2.3 stat ± 1.5 SJ/S fs. 

This measurement was averaged with previously published DELPHI results [2,4]: 

r r (l-prong, 91) = 298 ± 7 stat ± A sys fs, 
r r (l-prong, 92-93) = 291.8 ± 3.3 stot ± 2.1^ fs, 

to provide the result for all LEP-1 DELPHI data for the one-prong methods: 

r r (l-prong, 91-95) = 292.3 ± 1.8 stat ± l.2 sys fs. 

The final one-prong estimation was combined with the three-prong measurement to 
give the best estimation of the tau lifetime from the DELPHI data. Only the systematic 
uncertainty attributed to the alignment of the vertex detector is common between the 
one-prong and three-prong measurements, resulting in a 5% correlation between the two 
results. By combining the two results and taking into account this correlation, a tau 
lifetime of 

r r = 290.9 ±1.4 stat ± 1.0 SJ/S fs 

was obtained. 

5 Summary and Conclusions 

The tau lifetime has been measured using the DELPHI LEP-1 data sample. The result 

r T = 290.9 ± lA stat ± 1.0 sys fs 
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was obtained. This result supersedes all previous DELPHI measurements of the tau life- 
time. The measurement is compatible with the values published by other experiments [22] 
and has a slightly better precision. 

Tests of r — /j, and r — e universality can be performed using this result in conjunction 
with the published DELPHI values for the r leptonic branching fractions [23]: 

B(t~ - e-z/ e z/ T ) = (17.877±0.109 stat ±0.110 sys )%, 

B{ T - -> fi'i7^ T ) = (17.325 ± 0.095 sto ± 0.077 SJ/S )%, 

together with the world average values of the lepton masses and the muon lifetime [5]. 
Using equations (1) and (2), and accounting for small radiative corrections (see [23] for 
a discussion), yielded 

— = 1.0015 ±0.0053, 
il = 0.9997 ±0.0046. 

9e 

The branching fraction measurements contributed to the uncertainty in these estimates 
with ±0.0043 and ±0.0035 respectively. 

Under the assumption of e-/x universality, = g e = g C)fl , it was possible to give a 
more stringent test of universality of the coupling of the r and that of the two lighter 
leptons. The two measurements were combined into one leptonic branching fraction, B efM , 
correcting for the phase space suppression of B{r~ — > yTv^Vj): 

B e ^ = (17.838 ± 0.066 stat ± 0.068 SJ/S )%, 

to compare the tau charged current coupling to that of the two lighter leptons. The result 

— = 1.0004 ±0.0041 

9e,fi 

was obtained, in excellent agreement with r-(e,/i) universality. The relation between the 
leptonic branching ratio and the r lifetime is shown in Fig. 7, under the assumption of 
e-fj, universality. The band reflects the uncertainty on the tau mass [5]. 
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Figure 1: Decay length distribution for three-prong tau decays. The points are 1991-1995 
data sample and the histogram is the result of the fit. 
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Figure 2: Distribution of the difference of impact parameters versus the projected acopla- 
narity for the joint '94-'95 data sample (•) and simulation (o); the straight line is the 
best weighted x 2 fit on all the 94-95 events. 
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Figure 3: Upper panel: fitted decay length versus trim fraction: '94 simulation (•), '94 
data (o), '95 simulation (■) and '95 data (□). Lower panel: the difference between data 
and simulation for '94 (•) and '95 (■). Statistical uncertainties of the fit are added for 
comparison on two data points. The simulated decay length is L sim = 2237 /im. 
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Figure 4: Average value of the variance of the pull function of the reconstructed impact 
parameter for simulated r decays, as a function of the normalized momentum (see text) 
for hadrons (A), muons (■) and electrons (•). By construction the lifetime information 
is removed from this pull. 
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Figure 5: Distribution of the miss distance for the calibration samples of dimuon, dielec- 
tron and two-photon events from 1994 data (o) compared with simulation (•) and the 
parameterization used (histogram). 
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Figure 6: Distributions of the miss distance for the 1994 and 1995 data samples, points 
with error bars are data, the histogram the fitted distribution. The lower plots show 
the pulls of the differences between the observed data and the expected events from the 
simulation. 
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Figure 7: The measured values of the leptonic branching ratio and the lifetime of the tau 
lepton as measured by DELPHI (point), compared to the Standard Model relationship 
(shaded band). 



